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Abstract 



In this thesis, we introduce expansion systems as a general framework to 
describe a large variety of approximation algorithms, such as Taylor ap- 
proximation, decimal expansion and continued fraction. We consider some 
basic properties of expansion systems, and also study criteria for conver- 
gence. Further, we introduce the notion of isomorphisms between expansion 
systems. In the appendix, we discuss another class of expansion systems, 
which we call approximation systems. Many claims of convergence in this 
appendix, remain to be proven. 
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Introduction 



The concept of approximation is an important concept throughout the realm 
of mathematics. Many different approximation algorithms have been de- 
veloped, which are used in various contexts. But despite the variety of 
algorithms, there seem to be particular mechanisms which some of these 
algorithms have in common. For instance, many approximations are based 
on the existence of a fixed point for a certain contractive mapping, so that 
under the right conditions, iterated application of this mapping will approx- 
imate this fixed point. Much theory has been developed to capture this 
phenomenon, and its theorems are widely applicable (see for instance [13j). 

In this master thesis we will focus on a different approximation mecha- 
nism, which also seems to underlie many known approximation algorithms. 
It goes about as follows: we have an object y, on which we will apply a se- 
quence of mappings. After each of these mappings, our y has changed, and 
we record the essential information of this result before going to the next 
mapping (for example we repeatedly differentiate a function, and after each 
time we record its value in a point xq). Then, after having recorded these 
values up to a certain level, we try to reverse process all these mappings, 
using just the data we have collected so far (for example we repeatedly in- 
tegrate a function from the point xq, using the recorded values in xq as the 
constants of integration). In many examples, this algorithm turns out to 
give good approximations for the initial object y. The more of the recorded 
coefficients we involve in our approximation, the more accurate it generally 
becomes. This observation should motivate further investigation to this gen- 
eral principle, and we will refer to these algorithmic structures as expansion 
systems. 

Besides the example of Taylor series, which (as the reader may have fig- 
ured out) just served to illustrate the underlying mechanism of expansion 
systems, there are many more approximation algorithms that can be viewed 
in this perspective. There is the decimal expansion, Fourier series, contin- 
ued fraction, Engel expansion, Newton's forward difference formula, and so 
on. We will come across all of these examples as the theory of expansion 
systems is being unfolded. 

In Section 1 of the thesis, a formal description of expansion systems is 
given, which will serve as the main framework for all later sections. In Sec- 
tion 2, we will discuss some basic features of expansion systems regarding 
the approximations it provides. As we will see, the objects that arise as 
approximations will actually have much more meaning in the context of ex- 
pansion systems, than just being an approximation. In Section 3, we will 
develop the theory which is essential for investigating the convergence of 
expansion system. We will also prove some elementary convergence theo- 
rems, as well as a convergence result for the class of monotonia expansion 
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systems. The subject of Section 4 is isomorphisms between different expan- 
sion systems. The concept of isomorphisms may help us to translate results 
from one expansion system to another, whenever such an isomorphism can 
be established. 

In the final appendix, we will shed some light on a type of approxima- 
tions, called approximation systems. In a somewhat different setting, these 
type of approximations have been the subject of the author's bachelor the- 
sis Approximation Systems (in Dutch: Benaderingsstelsels) |18j . and they 
are also the subject of the identically named paper [19], coauthored by Prof. 
Koornwinder. In this thesis, we will reformulate this type of approximations 
in the context of expansion systems. As of yet, there is still a discrepancy 
between the results regarding approximation systems that have been estab- 
lish theoretically, and the results that have been experimentally obtained. 
For the treatment of the subject in the appendix, I have decided to espe- 
cially share the more practical results, accepting the lack of foundation to 
prove them. 

The subject of this thesis came into being during work on the article 
Approximation Systems. Prof. Koornwinder asked me to find out whether 
other types of approximation algorithms could somehow be fitted into the 
theory of approximation systems, and after much puzzling I decided that 
this was going to be difficult. However, there were nonetheless some basic 
principles which the approximation systems shared with other types of algo- 
rithms, and it became clear that there is in fact a general framework which 
all of these algorithms would fit in. When I decided to work out this idea as 
the subject of my master thesis, it was not quite clear to me yet, how much 
more I would be able to write down, besides the very definition and a couple 
of examples. But along the way, it appeared that there was actually quite a 
lot of theory to explore, and an initial habit of searching for content, soon 
enough changed into a habit of focusing on just a couple of basic features. I 
can therefore only apologize for all possible omissions, that probably would 
better have been included into this thesis. For example, like Prof. Homburg 
has remarked, it would make sense also to discuss the situation in which the 
expansions are linear maps. Although some of the examples in this thesis 
do belong to this category, such as the Taylor expansion and the Fourier 
series, a general theory for linear expansion systems is not given. And as 
Prof. Koornwinder suggested, also the theory of fast wavelet transforms 
could possibly be linked to expansion systems, but this hasn't made it into 
the final thesis either. Hopefully, future efforts will fill in gaps like these. 
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1 Description of Expansion Systems 



We consider the following situation: we have a sequence of so called element 
spaces S = {Si}i^fq, where N = 0, 1, . . ., and in each set Si there is a neutral 
element Vi ^ Si. For instance, if our Si are additive spaces, this will often be 
the zero element. Further we have another sequence of so called coefficient 
spaces C = {CjjjgN- Between these sets, we have the following collection of 
injections: Fi = {Pi,Ei) , where Pi: Si ^ Ci and Ei'. Si ^ 'S't+ij mapping 
the neutral element of Si to the neutral element of S'i+i, i.e. 



Ei{vi) = V, 



(1.1) 



Note that Pi and Ei themselves don't need to be injective. In diagram this 
looks like shown in Figure [T| 

-Eo El E2 Es 
JQ > Ol > D2 > 03 > • • • 





Pi 




P2 















P3 



Co 



Co 



C3 



Figure 1: Structure of expansion system 



Such a system 



{S,C,F) := {C'ilieN, {-fi}^ 



we call an expansion system, or ES for short j^] Each Pi in our ES will 
typically function as a projection on some simpler space, such as Z or C. 
The mapping Ei typically maps Si onto a similar space Sj+i, but thereby 
losing some information which is contained in the value of Pi. 

Example 1.1. For i E N let {Si,Vi) = ([0,1), 0), the set of real numbers 
greater or equal to and smaller than 1, with being the neutral element. 
Let = {0, 1, . . . , 9}. Furthermore, let Pj : 5j — >• Ci be the projection onto 
the first digit after the dot (i.e. Piy := [lOyJ, where [-J denotes the floor 
function), and let Ei: Si — )• Si+i be defined as Eiy := lOy — Piy. This 



expansion satisfies (1.1). Now Fi := {Pi,Ei) is injective, for the equations 
[lOyJ = [lOy'J , and lOy - [Wy\ = lOy' - [Wy' \ , 

Hence = (5, C, F) defines an expansion 



imply of course that y 
system. 



y 



*Whenever we refer to a sequence of element spaces S, we implicitly assume the exis- 
tence of a sequence i/i of neutral elements, even when not mentioned explicitly what these 
neutral elements are. 
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Note that in the example above, when we apply Ei to some element 
y E Si, we have lost some information about yi in the value Ei y, namely the 
value of its first digit after the dot. This information is then contained in 
the value Pj y. But despite our loss of information, Ei really is a surjective 
mapping from Si to iSj+i, which spaces are both equal to [0, 1). This is pos- 
sible because of the expanding nature of Ei. As this is the typical pattern 
in most ES's, wc will refer to the mappings Pj as our projections, and to the 
mappings Ei as our expansions. However, the mappings Pi and Ei do not 
necessarily have to obey this particular pattern, and the words projections 
and expansions are therefore merely terminology. 

Because all Pj are injective, we also have a sequence of inverses Ff^ : Fi{CiX 
Si) — >■ Si. Now given an y G Sq, we define yo := y and yj+i := P^yi, and 
furthermore Cj := Pj yi, so that we get yi = Fi^{yi^i, Ci). These q are called 
coefficients, and the full sequence c := (co,ci,...) is called the coefficient 
code. 

Remark 1.2. For convenience, we will adopt the special notation Pj for 
the accumulative composition of the first i expansions, so that we have 
Eiy = Ei ■ ■ ■ EQy = yi. Likewise, we write Pj for the composition of Pj with 
Pj, so that wc have Pi y = Pi Ei ■ ■ ■ Eq y = Ci. Also to the Pj : S'o — )■ Cj, we 
will refer to as projections, and the Cartesian product of all these projections 
gives rise to the so called coefficient map 

P:So^C, y^{co,ci,C2,...), 

where C := HieN ^« stands for the Cartesian product of all the coefficient 
spaces. In a certain sense, this map codifies elements of Sq as a list of 
coefficients from the spaces Cj, but it does so not necessarily injective. 

Having settled down the prerequisites of an ES it turns out that under 

the right conditions our ^ gives rise to an approximation for y G S, based 
on its coefficients c = {cq, ci, . . .) = P y . This approximation is constructed 
as follows: 

y|r' := (1.2) 
yl"] := Pj-^Cj, ygi) (i = n - 1, n - 2, . . . , 0) (1.3) 

provided that 

(Q,yJi) GPi(5j), for i = n-l,n-2,...,0. (1.4) 
Our ultimate n-th order approximation will then be: 

yW := yM (1.5) 

The above construction of y^'"'\ leads us to the following definitions: 
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Definition 1.3. An expansion system ^ is called proper at order n w.r.t. 
y £ Sq (or equivalently, to tiie coefficients cq, ci, . . . , c„_i), if it satisfies 



(1.4) tlirougliout the process of constructing y'-^K Only when dealing with 
an ES which is proper at order n, this construction will succeed, and are 
we allowed to talk about the n-th order approximation of y, also called the 
n-th convergent of y. When ^ is proper at any order n w.r.t. y, we say ^ is 
(completely) proper w.r.t. y. Also, when we say that ^ is proper without 
mentioning w.r.t. which element, it is assumed that it is proper w.r.t. any 
element in 5*0 . As this concept will return quite frequently, we will make use 
of special notation for the set of elements w.r.t. which ^ is proper: 

Sq"^ := {y G So ■ ^' is proper at order n w.r.t. y}. (1-6) 

Likewise, we let Sq~^^ and Sq"""^^ denote the elements in 5o, w.r.t. which 
^ is proper up to and including order n, respectively completely proper. 

In situations where we are dealing with multiple ES's at the same time, 
confusion might arise about relative to which ES we are determining the 
approximating In that case we will write yt*'"-! instead, to stress the 
ES we are working with. 

We will end this section with a simple proposition about properness, and 
give some applications to illustrate its use. 

Proposition 1.4. Suppose ^ is such that all Fi: Si ^ •S't+i x d are bijec- 
tive. Then ^ is proper. 

In] 

Proof. We have for z = n, n — 1, . . . , 1 that, given y\ £ Si and q^i G Ci_i, 
the surjectivity of implies (cj„i,y|"'') S Fj_i(S'i_i). Hence (1.4) 



IS 



always satisfied. □ 

Note that we could just as well have required Fi to be surjective, as we 
already have that it is injective by definition of an ES. This proposition 
implies for instance properness of the ES ^ in Example |1.1[ To see this, 
let's consider this ES once more here. 

Example 1.5 (Decimal Expansion). Let ^' be as in Example |1.1[ It 



is proper by Proposition 1.4, as we have that each Fi is bijective, and its 
inverse is given by: 

So if the coefficients of some y £ Sq are cq, ci, . . ., then we get as n-th order 
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approximation: 



10 

Cn-2 , Cn-1 



yn-2 — in 



n 



10 100 



Co Ci , Cn-l 



yo - 1 n + T nn + • • • + 



10 100 "■ 10"-l' 

In other words, the n-th order approximation coincides with the decimal 
number: 

yW = O.C0C1C2 . ..Cn-1 

It is good to reahze that it really requires a proof that this particular se- 
quence of decimal approximations is not only convergent, but that the num- 
ber it converges to is also the number we used as input of the algorithm, so 
that indeed we may suggestively write: 

y = O.C0C1C2C3 . . . 

Although an ad hoc proof of this wouldn't be that difficult to provide, we will 



postpone this to Theorem 3.22, where a more general result will be proven. 
Further it should be remarked that, although the space Sq is restricted to 
the set of non- negative numbers smaller than 1, we can quite easily adept 
the ES so that it can take bigger numbers as input. If we want Sq to be 
M>o, we just put Poy = [logig y] for non-zero y, and E^y = y/10^o(y) for 
non-zero y, and Fq{Q) = (0,0) for the case y = 0. This extends the ES to 
non-negative numbers of arbitrary size. 

Example 1.6 (Continued Fraction). Just as in Example |1.1[ we let 
{Si,i>i) = ([0,1), 0), but now we let d = Z>i U {00}. Furthermore, let 
Pi be given by: 

I 00 if y = 

I [1/yJ , otherwise, 



P^y 

and Ei be given by 

Eiy = 



if y = 

1/y — P-iy otherwise. 



This satisfies Equation |1.1[ and Fj = {Pi,Ei) is a bijection, for which we 
have as inverse: 

if c = 00 
otherwise. 



F-'ic,y) 
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By Proposition 1.4, ^ is proper and for its n-th approximation we find: 



1 



n 

Vn-l 



[n 
yn-2 



Cn-1 

1 



1 

Cn-2 + 



Cn-1 



Co + 



Cl + 



1 

... + — 



Cn-1 



It turns out that tliese rational approximations converge to the initial ele- 
ment y asn goes to infinity (cf. [Sj §1.6]), which justifies the concise notation: 

1 

y = co + (1.7) 

C1 + — 

C2 



C3 + . . . 

Just as for the decimal expansion, we will prove later on in this thesis that 
this is indeed valid for all y £ Sq. 
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2 Orders and Convergents 



In the previous section, we have seen how an expansion system codifies its 
elements in a sequence of coefficients. We have also discussed how these 
coefficients on its turn can give rise to an approximation. In this section, we 
will see that the approximations are characterized by the property that they 
are eventually mapped by the expansion sequence E to a neutral element 
I'n £ Sn- The values of its first n coefficients will coincide with the element 
it approximates. We will also see that the existence of elements with these 
properties, is at the same time a guarantee of properness. Throughout all 
the propositions in this section, we will assume we are given some expansion 
system and we again adopt the notation of the introduction for all of its 
components. 

Definition 2.1. We call y £ Sq an element of finite order, if for certain 
n G N we have yn = Vn- The smallest such n, we call the order of y. If 
y is not of finite order, we say it is of infinite order. As with the property 
of properness, we will use special notations Sg"^, ^q""^, •Sq^'^', and 5q°°' to 
denote the elements in Sq respectively of order n, order less than or equal 
n, finite order, and infinite order. 

Remark 2.2. By |1.1[ the order of an element y could also be characterized 
by saying: y has order n, iff y^ = z^m for all m > n. This is equivalent to 
the definition above. 

Theorem 2.3. We have y € s'^\ iff there is an element y' G S^q \ having 
the same coefficients up to and including order n — 1, i.e. Piy' = Piy = Ci 
for i <n — \. If this is the case, then y' is unique and equals yl"! |^ 

Proof. Proof of the forward implication: one verifies that it follows directly 
from the construction of y'"! that EiyY^ = and so we have for the 

composition Ei of expansions: 

(yM), = i,yM=i^^...^„yN=yM (2.1) 

and thus in particular (yf"^)n = yn"' = Vn, which means y^"! G S^q \ More- 



over, (2.1) gives us that: 



A = P^ {y^% = P^ yf ^ = P^ F{\yf},,ci) = c, = h y. 

Hence yt"'' satisfies the requirements for our y'. 

Now we prove the backward implication. We just carry out the construction 
of y["', and let ylT^ = Because we know that y^_i is an element in Sn-i 



^Throughout this thesis, primes are only used to distinguish symbols (e.g. y and y') 
and do never have the meaning of differentiation. 
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for which Fn-iy'n-i = (^'n,c„_i), it follows that = Ff^ {'i}n\ Cn-i) 

is well defined, and equals y'^-i- The same observation now applies to 

Fn-iy'n-2 = (yn-i,c„-2), SO that y^^l^ = y'n-2- Continuing this process, 
we get at the last step that t/q"' = yt"' is well defined and equals y'^ = y'. 
Because of the arbitrarity of y', and the fact that ^q"' is uniquely determined 
by its construction, the additional remark follows that this y' is unique. □ 

Corollary 2.4. Given y € we have that y € 5'q~"\ and y'™'' = y for 
all m > n. 



Proof. This follows from the observation that by Remark 2.2 = i^m for 
all m > n. So for each order m > n, y itself serves as the looked after 
element z, and so we have yt*"] = y. □ 



Remark 2.5. Theorem 2.3 not only provides another formulation for proper- 
ness, it also characterizes the n-th convergent of an element y as the 
unique n-th order element z in 5o, having the same coefficients as y up to 
and including order n — 1. The corollary on its turn, implies that 

4"' = {^'"^ ■■ y e 4"^}- 

In particular we have that si"^ C 5q~"\ 



We can apply Theorem 2.3 to the following example of an expansion 
system: 

Example 2.6. Taylor Series 

Let each Si be the space of holomorphic functions on some open neighbor- 
hood [/ C C of a point xq G C, and let Cj = C and z/j the zero function. 
Furthermore, given a function y{x) G Si, let Piy = y{xo) and 

TP ( X y{^) - Piy 

Eiy{x) = , 

X — Xq 

the value in xq being defined by its limit. One verifies that for given y{x) = 
Co + ci{x — xq) + . . ., we have Piy = Ci. We have that F, = {Pi, Ei) satisfies 
and furthermore ^ = {S, C, F) is proper, for given an order n, we have 
that z = Co + ci{x — Xq) + . . . + Cn-i{x — xo)""^ satisfies Theorem 



2.3 



(We could also have invoked Proposition 1.4, noting that every holomorphic 
function has an integral with prescribed value in the point xq.) In this ES, 
the elements of order n are just the polynomials of order n — 1 (the zero 
polynomial in case n = 0) . The elements of infinite order are exactly all non- 
polynomial holomorphic functions. The n-th order approximation coincides 
with the n — 1-th Taylor approximation at the point xq. 
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It is important to realize that if an ES ^ is proper w.r.t. y at order n, this 
doesn't necessarily mean that ^ is also proper w.r.t. y at orders below n, 
nor that it is proper at orders above n. We will use a variation on Example 



2.6 as a counterexample. 



Example 2.7. Again we consider the expansion system ^ as given in Ex- 
ample 2.6, but this time we let xq = and U be the closed interval [0, 1]. 



Further, we let Vi,Ci,Pi and Ei just be as in Example |2.6[ but we restrict 
each Si only to those elements yj, such that \\yj\\ < 1 for all j > i, i.e. 
yi remains in absolute value < 1 after composition with arbitrarily many 
expansion Ej, j = i,i + 1, . . .. In particular, the elements in Si itself are of 
absolute value < 1. Note that by this condition, we still have Ei{Si) C 5'i+i. 
Now consider the polynomial y = ^ + x — x"^ + — x"^. We can calculate 
that: 



yo 

yi 
y2 
ys 
y4 
y5 



+ X 

- X + X 

1 + x- 

- x\\ 



x^ + x^ 



X 



0.60583 < 1 



1 



ye I 



loll = 0. 



All these yi have absolute value < 1 on the interval [0, 1], and so yi G Si. It 
turns out that y S ^q^^ , as we have ior z = \ ^ x — x^: 



Fol 
ksl 



,1 2, 
|- + X - X 

|1 -x|| = 1 

I - lll = 1 



< 1 



loll = 0. 



So z is an element in and it satisfies Theorem 



2.3 



by which we may 

conclude properness at order 3 w.r.t. y. But for order n = 2, our only 
+ X, but II i + x|| = li > 1, so there can be no such z in 



candidate is z 



So- Likewise, for order n 
but 



\l+x 



x^ + x^l 



4, we have 2 + ^ 



x^ + x^ as only candidate, 



1^ > 1, so this z does not lie in 5o either. Hence 



^ is proper w.r.t. y at order 3, but not at order 2 and 4. 



(n) 

Proposition 2.8. Suppose y G •S'g , 
coefficients as y, i.e. Piy = Piy' for < i 



and we have yl""' = y' 



An] 



and y' G Sq has the same first n 
< n — 1. Then also y' S 'S'q"^ 



Proof. By Theorem 2.3 we conclude that there is an element y^"' E S^ with 



the same first n coefficients as y. So by assumption, yt"J also has the same 
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first n coefficients as y' . Therefore, again by Theorem 2.3 we have that 

□ 



y' G 5j"\ and y^" 



y 



'[n] 



Corollary 2.9. Suppose for certain n and m < n we have y ^ Sn H S, 



(n) <^(m) 



I.e. 



is proper w.r.t. y both at order n andm. Then also y^^' S Sq n5, 



^0 

(n) n{m) 



and we have {y^^' 



■"0 



I.e. 



the m-th order convergent of y'"' equals 



the m-th order convergent of y itself. In particular (yN)N = yb 



Proof. By Theorem 2.3 , y and yl"I have equal coefficients up to and including 
order n — 1. Then certainly y and yt"' have equal coefficients up to and 
including order m — 1. Hence Proposition 2.8 applies, which implies that 

y[n] g Q{m)^ (yN)M ^ yH. □ 
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3 Convergence of Expansion Systems 



When we are given an ES ^ with a topological space 5o, it would be nice 
if y["l becomes a better approximation of y as n increases, so that we can 
define 2/'°°^ := lim y'""' = y. However, this may not always be the case, and 
the following definition provides names for the situations that may arise: 

Definition 3.1. Suppose ^ is proper with respect to y in the topological 
space 5*0 . We then call ^ convergent with respect to y, if y^"''^ converges 
to some element in So, for which we will use the notation y^°°h We say 
^ is accurate w.r.t. y, if it is not only convergent with some limit y^°°\ 
but moreover equals y. We say ^' is divergent with respect to y, if 
the sequence of convergents y'"' does not converge at all. Furthermore, we 
simply say ^ is convergent or accurate, if it has this respective property 
with respect to all of its elements. 

Remark 3.2. Note that by Corollary |2.4[ an ES is always accurate w.r.t. its 
elements of finite order, even if Sq is given the discrete topology. 

In this section we will investigate some of the conditions under which ^ 
is (accurately) convergent. As convergence depends on the topology which 
the element spaces are endowed with, it is important to pay attention to 
the topological aspects of ES's. As it will turn out, in many cases the 
approximation will converge in one topology, but not in the other. We will 
therefore first develop a framework which will help us in investigating these 
phenomena. 

3.1 Topologies on Expansion Systems 

Definition 3.3. Let f : X ^ Y he a map between X and Y, where Y is 
given the topology '3^. Then we define the generated topology /'^(^), to be 
the topology on X generated by the map /, meaning it is generated by the 
sets f-'^{U), where U is open in ^ (cf. [HI 1.4.8]). 

More specific, let 'if be a topology on the space C of coefficient sequences. 
We then define P~^{'rf) to be the topology on Sq, generated by the sets 
P~^{U) with U open in and call this the generated topology. Usually 
each space Cj already has its own topology ^i, and unless otherwise spec- 
ified, we will assume that is the product topology with respect to these 
topologies. In that case, we call the generated topology P"^(^) the weak 
topology generated by the collection of projections Pi to the spaces {Ci,%) 
(cf. [n, 3.16]). In the context of expansion systems, the weak topology on 
^0 which is generated by the discrete topologies on Ci turns out to be an 
important topology. We will call this topology on 5*0 the weakly discrete 
topology, and denote it by P"^(i^*), where S>* on its turn denotes the dis- 
crete product topology on C. The discrete topology on C will be written as 
just 2. 
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Remark 3.4. The generated topology P~^{^) can also be defined as the 
weakest topology such that P: 5o — )■ C* is continuous. In particular, the 
weak topology is the weakest topology such that all Pj : 5o — )• are contin- 
uous (this follows from the fact that P is continuous iff each Pi is continu- 
ous, see [I6l 19.6]). This topology is therefore generated by sets of the form 
P-^{Ui), where i £ N, and Ui is open in Cj. This implies that a sequence 
yH in 5*0 converges to y in the weak topology, iff for each i G N: Pi 
converges to Pm in the topology 'rfi of Cj. Equivalently, yl*^] converges to 
y in the weak topology, iff PyW converges to in the product topology 
^. The weakly discrete topology P~^{^*) is likewise generated by the sets 
Pf^{ci), where i £ N and Cj runs through all the points in Cj. In particular, 
yH converges to y in P~^{&*), iff for each i S N: Piy^"^^ equals Piy for n 
large enough. The weakly discrete topology is the strongest one among all 
possible weak topologies generated by the Pi, so given an arbitrary weak 
topology, the weakly discrete topology is necessarily stronger or equal to it. 
Although it is useful to use different names for the generated topology and 
the weak topology, we should keep in mind that in a strict sense there is no 
essential distinction between these two types of topology, because a weak 
topology can be viewed as the topology generated by the product map to 
the product topology, and a generated topology can be viewed as a weak 
topology for just one map. 

Lemma 3.5. Consider an ES ^ where Bi {i G N) are bases for the topologies 
'^i of Ci. Then the product topology '^o on C has a basis of the form: 

TT:^{Bo, Bn-l) = BqX ...X Bn-l X Cn X Cn+1 X . . . , (3.1) 

where Bi G Bi and TTn : C* — )• Yli^n ^« projection on the first n coordi- 

nates. The weak topology P~'^{'i§') on Sq has as a basis: 

n-l 

p-\Bo,...,Bn^i):= f]P-\Bi), (3.2) 

2 = 

where Bi £ Bi. In particular, we have that Ss* has a basis of the form: 

7r„^(co, . . . , C„_i) = {co} X ... X {Cn-l} xCnX Cn+1 X . . . , (3.3) 

and P~^{&*) has a basis: 

n-l 

p-\co,...,Cn-i)= f]Pi\c,), (3.4) 
i=0 

where each Ci is a point in Ci. 

Proof. We recall that the product topology 'rf is generated by all sets of the 
form n ieN ^« ' where for only finitely many i G N we have Cj / S Bi (see 
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for instancefTF, Thm. 19.2][5l 2.3.1]). Given such a set HieN^' ^® 
assume that 14 = Cj for z > n. Now let c = (cq, ci, . . .) G OieN each 
z < n we pick a basis element Bi of ^Bj, such that Ci £ Bi C Vi (when Vi = Ci, 
any Bi containing q will suffice, and when Vi ^ Ci,we can just set Bi = Vi). 
Now 7r„^(i?o, • • • , Bn-i) is an open set containing c, and which lies within 



the set HieN ■'^^ arbitrarity of HieN ^'^^ ^' this proves that (3.1 ) 



defines indeed a basis for Equation (3.2) follows from the observation 



that if ;S is a basis of then {P~^{B) : B £ B} serves as a basis for 



Equation (3.3) and (3.4) are now special cases of respectively (3.1) 
and ( |3.2[ ), where we note that the collection of all singletons forms a basis 
for the discrete topology. □ 

3.2 Some Elementary Convergence Theorems 

Theorem 3.6. Given y G Sq'"°°\ we have that Py^"^^ converges to Py in 
the discrete product topology Qi* . Further more, y^"^ converges to the initial 
element y £ Sq in the weakly discrete topology 



3.5 



Proof. Let y be a neighborhood of P y in the topology . By Lemma 
without loss of generality we may assume that V is an element of the form 



7r~^{cQ, . . . ,Cn-i). By Remark 2.5 we have that Piy'^^' = Piy = Ci for all 
i < m, and so in particular for all i < n, when m > n. Hence g V for 

m > n, which proves by arbitrarity of V that Pyl™-! converges to Py. The 



second claim follows by Remark 3.4 □ 



Convergence in the weakly discrete topology turns out to be useful for 
our understanding of expansion systems, but its direct implications for the 
convergence of yl"^ in other topologies are nonetheless limited. The following 
example will demonstrate this. 

Example 3.7 (Fourier Series). Let Si be the space of complex integrable 
functions on [0,27r], such that /(O) = /(27r). For Ci we take the set C^. 
Now we define 

Py:={ r e'''y{t)dt, r e-'''y{t)dt) = {c-„c+i), 
Jo Jo 

where i := Furthermore, we set 

Ei y{x) := y{x) - C-iC'"''' - c+iC 



ILX 



One verifies that Piy = (c_i,c+j) corresponds with the — i-th and +i-ih 
Fourier coefficient of y. As the n — 1-th order Fourier approximation y'^ := 
c_(„„i)e~("~^)^^' + . . . + c„_ie("~^)^^ is clearly an element of 5o, and we have 



that En{y'n) = = M„, it follows by Theorem 2.3 that y^*^' equals the n — 1-th 
order Fourier approximation. 
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Here, Theorem 3.6 does not guarantee that yl^J will converge in other 
topologies than the weakly discrete topology. Indeed, in the example of 
Fourier approximation, a function may have well-defined Fourier coefficients, 
so that y'"' converges in the weakly discrete topology, while its Fourier series 
diverges in the topology of pointwise convergence (cf. [22l Ch. 3, §2.2]). 
In fact, the weakly discrete topology is in this ES not even weaker than 
pointwise convergence, but just incomparable. For if we take the sequence 

/„:[0,2^]^C, fn{x) = l/n, 

it clearly converges pointwise, but its zeroth Fourier coefficient does not 
converge in the discrete topology. However, given a topology we have in 
mind for 5o, it is often possible to turn the ES into a convergent one, simply 
by restricting the space 5*0 to a smaller space. For example, in the case of 
Fourier approximation, in order to get pointwise convergence it suffices to 
restrict 5o to differentiable functions (see [22, Ch. 3, §2.1]). 



Another important thing to notice, is the following. Theorem 3.6 tells 
us that converges to y in the weakly discrete topology. But, when we 
alter the values of the function y on a set with Lebesgue measure 0, then 
this will not have any effect on the Fourier coefficients Cj, and therefore 
also the convergents y'"' will not get altered. But according to Theorem 



3.6, should also be accurately convergent towards the altered y in the 
weakly discrete topology. This means that P~^{&*) is not Hausdorff, and 
yM converges to multiple different elements at the same time. It turns out 
that the property described in the following definition, is characteristic for 
the weakly discrete topology to be Hausdorff. 

Definition 3.8. An ES ^ is called separating, if the map P: 5*0 — )• C* is 
injective. In other words, ^ is called separating, if the family of projections 
Pi separates points in 5o, i.e. for y, y' G Sq and y ^ y', we have Piy ^ Pi y' 
for some i £N. 

So if ^ is a separating ES, the elements in can be uniquely codified 
as a list of coefficients (cq, ci, . . .) G C. 

Proposition 3.9. Let 'if be a Hausdorff topology on C , and let the ele- 
ment space Sq be endowed with the generated topology P~^{'if). Then Sq is 
Hausdorff iff ^ is separating. 

Proof. If ^ is separating, then the map P: Sq ^ C is a injective continuous 
map into a Hausdorff space, which implies that is Hausdorff. On the 
other hand, if ^ is not separating, then we have c = Py = Py' G C for 
certain y, y' G ^o, y ^ y'. Now if Vy and V^/ are two neighborhoods of y and 
y' in 5*0, we may w.l.o.g. assume that Vy = P~^{Uy) and Vyt = P~^{Uyi) 
for some open neighborhoods Uy and Uyi in C. Then c £ Uy Ci Uy', and 
so P"^(c) G Vy n Vy/, which shows that any two neighborhoods of y and y' 
cannot be disjoint. □ 
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Theorem 3.10. Suppose is a separating ES, such that Sq is a topological 
space and all Ci are topological Hausdorff spaces. Moreover, suppose ^ is 
convergent w.r.t y €z Sq. If P is continuous at y[°°] in the weak topology, 
then y = y^°°^ , which means ^ is accurate w.r.t. y. In particular, if ^ is a 
convergent ES and P is continuous on Sq, then ^ is accurate. 

Proof. If P is continuous at we have that Py^"^^ converges to Py^°°^, 



and thus by Remark 3.4, y'"' converges to in the weak topology. On the 



other hand y'"' converges to y in the weakly discrete topology by Theorem 



3.6, and thus also in the weak topology. As Sq is Hausdorff in the weak 



topology by Proposition 3.9, we must have that y = y^°^^. □ 



Remark 3.11. Note that in the above theorem, for P to be continuous at an 
element y^°°\ it is sufficient if all the expansions -Ej : — )• /Sj+i and all the 
projections Pi: Si ^ Ci are continuous at , as each Pi is composed from 
these maps. Also note that the theorem does not state that y'"' converges 
to y. It only says that if the sequence converges, and P is continuous at its 
limit, then this limit can only be y. 

Example 3.12 (Convergence of Fourier Series). We demonstrate how 
Theorem 3.10| can be applied to obtain accurateness of particular Fourier 



approximations. Consider the ES ^' from Example 3.7, except that we let 
the Si consist of only continuous functions, with the topology of uniform 
convergence. The projections Pi are then continuous on Sq, as we have the 
estimation: 

Piy-Piy' = Pi{y - y') = / {y{t) - y'{t))e'''dt 

Jo 

< / \yit)-y'{t)\dt<27r\\y-y'\\. 
Jo 

Using the property that Fourier coefficients are unique for continuous func- 
tions (cf. [22l Cor. 2.2]), we know that ^' is also separating. Now suppose 
y £ So has absolutely convergent Fourier coefficients q, i.e. 



E 



Cj < oo. 



Then one easily verifies that the series: 

n-l 



yN(^)= ^ 



CkC 



converges uniformly on [0, 27r]. Hence by Theorem 3.10 the approximations 
will converge to the initial element y. 
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Lemma 3.13. Let f : X ^ Y be a mapping, and let J^T be the topology on 
X. Then the following statements are equivalent: 

I ^ is weaker or equal to the generated topology f~^{&), where '2 is the 
discrete topology on Y . 

II Every open set U of ^ can be written as f~^{V), where V is a subset 
ofY. 

Ill f{p) = f{p') implies thatp andp' are topologically indistinguishable in 
{X, i.e. every open set which contains one of these two elements, 
also contains the other. 

Proof. The equivalency between I and II is straightforward. So let's consider 
II implies III. Then, if f[p) = f{p'), and we have a open neighborhood U 
containing p, it follows that U = f~^{V) for some V C Y. But then f{p) G V, 
so that p' £ f'^ifip)) C f~^{V) = U. By symmetry it follows that p and p' 
are indistinguishable. 

Now we prove III implies II. Now let U be an open set in X, and consider 
V := f{U). If p G f~^{y): than f{p) = f{p') for some p' G U , and so p and 
p' are topologically indistinguishable. Hence, we must have that p G C/ as 
well, which proves that U C f~^{V) C U. So U can indeed be written as 
f-^{V),withVcY. □ 

This lemma, together with the following definition, allows us to formulate 
another convergence theorem. 

Definition 3.14. Let f : X ^ Y he a map between topological spaces. We 
call such a map closed at a point y £Y, ii for every open neighborhood U of 
the set f~^{y), there is an open neighborhood V of y, such that f'^iV) C U 
(cf. P 4.5.13]). 

Remark 3.15. One can prove that a function f : X ^ Y is closed, iff / is 
closed at any point in Y according to the definition above (see for instance 

m 1.4.13]). 

Theorem 3.16. Let ^ be a proper ES, and let C be endowed with the 
discrete product topology . Assume So has a topology which is weaker or 
equal to P~^{&). Then ^ is acurate w.r.t. y, if P is closed at the point 
c = P{y) G C. 



Proof. Let f7 be a neighborhood of y. By Lemma 3.13 we have that U is also 



a neighborhood of P"^(c). So by assui nptio n, there is an open set V C C, 



such that P~^{V) C U, and by Lemma 3.5 we may assume this V is of the 
form 7r„^(co, . . . , Cn-i), where Cj = Pi y. As Pi y^"^^ = Ci for i < m, we have 
that yt"*! G P~^{V) C U for all m> n, which proves convergence to y. □ 
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Corollary 3.17. Let ^ be a proper ES, and let C he endowed with the 
discrete product topology Si* . Then ^ is accurate, if P is a closed map. 



Proof. This follows immediately from Remark 3.15 □ 



3.3 Convergence of Monotonia Expansion Systems 

In this subsection, we will define the so called monotonic expansion systems. 
This generalizes many instances of expansion systems, in which the element 
spaces are linearly ordered sets. The results we obtain for monotonic ES's 
can thus be applied to all the different instances of monotonic ES's. 

Before we come to the definition of monotonic expansion systems, we 
will first briefly go through some definitions related to orders on a set. We 
call a function f : X ^ Y between linearly ordered sets (X, >) and (y, >) 
monotonic, if / is either strictly increasing or strictly decreasing]^ Fur- 
ther, with the dictionary order on a (possibly infinite) Cartesian product 
X = X Xi X X2 X . . ., we mean the order such that (xq, xi, X2, . . .) < 
(x'o, x'^, X2, . . .), iff at the first entry k on which the two sequences differ, we 
have that < x^ (cf. Ch.l,§3]). The order topology on an linearly 
ordered space X, can be defined as the topology generated by the subbasis 
of open rays X^x' and X^x", with x',x" E X. In particular this topology 
has as a basis the collection of open intervals Xx'<-<x" with x',x" G X, 
together with the sets Xx<.<:x" and Xx'<:.<x in case X contains a minimum 
element x or a maximum element x (cf. [IS Ch.2,§14],[r, Ch. 1, §2.1])|§] A 
subset U of an ordered set X is said to be densely ordered, if for any two 
different elements x,x' £ X there lies an element u £ U in between. 

Definition 3.18. Suppose we have an expansion system ^, for which the 
sequences S and C both consist of only linearly ordered sets. Then, if each 
Fi : Si ^ CiX is monotonic when Cj x Si+i is given the dictionary order, 
we call ^ a monotonic expansion system. More concrete, this condition on 
Fi means that either for all elements y < y' in Sj, we have that: 

Piy<Piy' and (3.5) 
Eiy < Ei y whenever Piy = Piy . (3.6) 

or for all elements y < y' we have that: 

Piy>Piy' and (3.7) 
Eiy > Ei y' whenever Piy = Pi y' . (3-8) 



If all Fi are strictly increasing (i.e. we have (3.5) and (3.6)), we call ^ a 
strictly increasing expansion system. 



"^This meaning of the word monotonic might difler from those in other sources, where 
it could merely stand for strictly increasing or non-decreasing. 

^To be fully correct, we should remark that a basis must also contain the set X^<.<^, 
in the rare case that X only contains the single point x = x. 
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Example 3.19. The decimal expansion (Example 1.5), is a strictly increas- 



ing ES. One easily verifies that (3.5) and i3.6v both hold. The continued 



fraction (Example 1.6) is not a strictly increasing ES, but it is a monotonic 



ES, for we have that (3.7) and (3.8) holds for all « G N, where we should 



remark that oo is to be interpreted as the largest element in Cj. 

From the definition of a strictly increasing ES, we almost immediately 
obtain the following property: 

Proposition 3.20. Let ^ be strictly increasing expansion system, and let 
C be given the dictionary order. Then P : Sq ^ C is a non- decreasing map. 

Proof. Suppose we have Py > P y' for certain y < y' . Then there must be a 
k > 0, such that Piy = Pi y' for i < k, and Pf^y > Pkv' ■ But this contradicts 



(3.5) for i = k, from which we conclude that P is non-decreasing. □ 

This simple proposition will be useful in the proof of the following two 
important theorems. 

Theorem 3.21. Suppose ^ is a monotonic expansion system. If the set of 
elements of finite order 3^^°°'^ is densely ordered in Sq, then ^ is separating. 

Proof. First we prove the theorem for strictly increasing expansion systems, 
and then we will show how the result can be naturally extended to monotonic 
ones. So let us be given a strictly increasing ES ^, and give C the dictionary 
order. If ^ would not be separating, then there a re two different elements 
y'jV" € 5*0, such that Py' = Py". By Proposition 



3.20 



we know that P is 

non-decreasing, and hence Py = Py' for all y' < y < y" ■ In particular, we 
can pick two different elements z,z' of finite order, lying both between y' and 
y" (we can first apply the property of being densely ordered on y' and y" to 
get z, and then on z and y" to get z'). But then z and z' must have exactly 
the same coefficients, which contradicts that an element of finite order is 



uniquely determined by its coefficients (see Remark 2.5). Hence ^' must be 
separating. 

Now to show that the same is true if all Pi are only monotonic. This 
can be achieved by transforming our monotonic ^ into a strictly increasing 
expansion system ^' , by reversing the order relations on some of the spaces 
Si and Cj. Note that reversing an order on a linearly ordered space keeps 
the space linearly ordered. The reversing procedure will be according to the 
following rule: starting from z = 0, 1, . . ., we reverse the order on Si+i and 
Cj+i either if Fi is a strictly decreasing map in ^ and we haven't reversed 
the order of Si and Ci in or if Fi is strictly increasing in ^, but we also 
reversed the order on Si and Ci in ^' . By doing this, we keep al strictly 
increasing maps strictly increasing, and at the same time strictly decreasing 
maps in ^ become strictly increasing ones in ^' , so that our new ^' is a 
strictly increasing ES. As reversing orderings has neither any effect on the 
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element and coefficient spaces themselves, nor on the collection of expansions 
and projections (in fact we are only changing properties which lie outside the 
essential structure of the ES), the claim of the theorem also applies to ^''j^ 
Hence, the statement extends without problems to monotonic expansion 
systems. □ 

Theorem 3.22. Suppose ^ is a proper monotonic expansion system, and 
So has the order topology. Then ^ is accurate, iff ^ is separating. 



Proof. As in Theorem 3.21, we first prove the statement for strictly increas- 
ing ^, and then extend it to the case ^ is monotonic. 

First let us assume that ^ is accurate, but not separating. Then there 
are y,y' E Sq, having equal coefficients and is accurate w.r.t. both of 
them. This means that y'"! = y't"! converges to two elements at the same 
time, and so So cannot be a Hausdorff space. But 5*0 has the order topology, 
which is clearly an Hausdorff space. Hence ^ must be separating. 

To prove the reverse implication, ffist we give C the dictionary order, 
as well as the order topology ^ generated by this ordering. We now prove 
that the discrete product topology , is stronger than the dictionary order 
topology 'tf on C. For let V be an open set in the dictionary order topology, 
containing an element c = (co,ci, . . .). W.l.o.g. we may assume that V is 
either Cc!<-<c" or Cc<.<^~c" or C^^.^ji-, where c and c denote respectively the 
smallest and the largest element in C, if it exists. First we regard the case 
C'c'< <c"- As c is supposed to lie between c' and c", we have on the one hand 
that c' and c have equal coefficients up to and including n' — 1 and c^, < c„/ . 
On the other hand we have that c and c" have equal coefficients up to and 
including n" — 1, and c"„ < c„//. This means that any coefficient sequence 
with the same coefficients as c up to and including m = max(n',n") will 
lie in the set (?£'<■<£"• Hence vr^^(co, . . . ,Cm) is a neighborhood of c in the 
discrete product topology, that is contained in V . 

If y is a set of the form (?£<.<£", then we may assume that c is the 
minimal element c, for otherwise we can apply the previous result to 
But in that case, we can consider the set ■k~^{co, ■ ■ ■ ,Cm), where m is the 
fo'st coefficient such that Cm < c'^- A similar argument applies to the 
case V = Cgi^.<^-^, and thus it follows that the discrete product topology is 
stronger or equa l to the dictionary order topology. 



By Theorem 3.6, we have that Py^"J converges to Py in the discrete 
product topology. Therefore, PyN also converges to Py in the dictionary 
order topology This, in turn, implies thatw'""' converges to y in the 
generated topology P"^('^). By Proposition 



we know that P is non- 
decreasing, and as ^ is supposed to be separating, so that P is injective. 



3.20 



^Also the property of 5''^°°' being densely ordered in So cannot have been affected by 
the process of reversing orderings, even if we would have reversed the order on So itself, 
which we haven't. 
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we conclude that it must be strictly increasing. Therefore, we have that the 
topology P~^{^) is stronger or equal to the order topology on 5*0, for we can 
write (5'o)>y as P~^{C-^p^y,^), and likewise for (<S'o)<y". Hence, we have at 

last that yt"! converges to y in the order topology of 5*0, which means that 
^ is accurate. 

To show that the same is true if all Pi are only monotonic, we apply 



exactly the same trick as we have done in the proof of Theorem 3.21[ Again, 
this doesn't affect the essential properties of nor does it affect the order 
topology on Sq. Hence, the result follows for proper monotonic ES's as 
well. □ 

Corollary 3.23. Suppose ^ is a proper monotonic expansion system, and 
Sq has the order topology. If S^^°°^ is densely ordered in Sq, then ^ is 
accurate. 



Proof. This follows immediately by combining Theorem 3.21 with Theorem 



[3:221 □ 

This corollary can be used to prove convergence of the ES's in Example 

EH 



Example 3.24. The decimal expansion (Example 1.5) and the continued 



fraction (Example 1.6) are accurate, if Sq is given the standard topology 



(i.e. the Euclidean topology) on the real numbers. 

Proof. That these expansion systems are proper and monotonic, we have 
already shown in former examples. Noting that the standard topology on 
the reals coincides with the order topology (cf. \16\ §14, Ex.1]), it remains 
to prove that in both cases is densely ordered in Sq. For the decimal 

expansion, this is easy, as we have that the elements of finite order are exactly 
all multiples of lO"'^ between (and including) and 1, where k can be any 
positive number. These numbers are easily seen to be densely ordered (given 
two of them, we can multiply by a suitable power of 10, so that it boils down 
to picking a natural number between two other natural numbers). Hence ^ 
is accurate in the case of decimal expansion. 

In case of the continued fractions, we claim that the elements of finite 
order are exactly the fractions between (and including) and 1. To see this, 
we can look at the function M : Qo<-<i N, which assigns to a fraction 
p/q (where p and q are written in their lowest terms) the value p + q. Then 
one verifies that for {p,q) 7^ (0, 1): 

M(i5;,(^)) =m(^-^) =q-ic,-l)p<p + q = M(^^^ 

So if p/q is non-zero, then the value of M decreases with each application 
of Ei, while it should remain positive. This can only be the case if for 
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each fraction y in 5o, En y equals = j for n big enough, and this n is by 
definition the order of y. Also the rational numbers are of course densely 
ordered, and so by Corollary |3. 23 we have that the ES of continued fractions 
is accurate. □ 

There are many more instances of monotonic expansion systems. We will 
give some more examples. 

Example 3.25 (Egyptian Fractions and Engel Expansion). Let (5j, Vi) 
([0, 1),0), and Ci = Z>i U {oo}. Furthermore, let Pi be given by 



and Ei by 



P^y 



Ery 



oo if y = 

otherwise, 



if y = 

y — -p— otherwise. 



Now Fi = {Pi,Ei) is a bijection, and we have as inverse: 



F^\c,y) 



+ y 



if c = oo 
otherwise. 



By Proposition 1.4 this ^ is proper. If we give each Ci = Z>i U {oo} the 



reversed order, then each Fi is a strictly increasing map. In the way we have 
demonstrated in Example 3.24, we can prove that S'q'"°°^ is again exactly the 
set of rational numbers in Sq (this time, look at M{p/q) = p), so that ^ is 



accurate by Corollary 3.23, One easily sees that the approximation this ES 
provides, is just the sum of reciprocals of the coefficients. For example, if 
we take the element y = l/\/2 E 5*0, then we get the approximation: 



1 _ 1 1 1 1 

71~2"^5"^14l"^ 68575 



+ .... 



Breaking off the series at a given depth, the resulting approximation will be 
an expression which is a sum of only unit fractions. Such expressions are 
usually referred to as Egyptian fractions. 

Now let's make a slight change to our ^': we keep all components the 
same, except for Ei, for which we set: 



Eiy 







if y = 



y ■ Pi{y) — 1 otherwise. 
The inverse of Fi is then given by: 

if c = oo 



Fiic,y) 



l+y 
c 



otherwise. 
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Again, ^ is proper by Proposition 1.4, and each Fi is strictly increasing if we 
give Ci the reversed order. 51^°°' is the set of rationals (look at M{p/q) = p), 
and thus ^ is accurate by Corollary 3.23 For the n-th approximation we 
find: 



[n] 
Vn-l 





1 

Cn-1 



yn-2 



1 + 



Cn-1 



Cn-2 



1 + 

1 + 

1 + 



1 



Cn-1 



N ci 

% = 

Co 

Hence we have the concise expression for y: 

1 + 



1 + ^ 



ci 111 , , 

y = ^ = - + + + . . . (3.9) 

Co Co CoCi C0C1C2 

Here, the approximations are a special kind of Egyptian fractions, in which 
the denominators of consecutive fractions increase by a non-decreasing se- 
quence of factors cq, ci, . . .. Such approximations are called Engel expansions 

(cf. [Tirsj). 



Example 3.26 (Taylor Series). Let ^ be as in Example 2.6 except that 
we confine the spaces Si to only real holomorphic functions. We can view 
^ now as a strictly increasing ES, if we give Si the ordering defined by: 

y < y' if y{x) < y'{x) for all x £ {xo,xo + e), (3.10) 

for some sufficiently small e > 0, depending on y and y' . By the limit 
properties of y in xq, this ordering coincides with the dictionary order on the 
derivatives of y, and thus the Taylor coefficients. So this demonstrates that 
^ is monotonic. As for the accurateness of this is not a very interesting 
question, as Taylor approximations will at least locally converge just by the 
definition of holomorphic functions. But indeed, the elements of finite order 
are the polynomials, and in the ordering we have just put on our element 
spaces, these are densely ordered in the set of holomorphic functions. So 
by Corollary |3.23[ the Taylor polynomials converge to the initial function in 



the ordering given by (3.10). 
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The monotonic expansion systems seem to capture quite some well-known 
expansions. As for the continued fraction, in which the expanding maps Ei 
are all strictly decreasing, a similar generalization has priorly been studied 
by Bissinger in [3j, and for the decimal expansions, in which the expanding 
maps are strictly increasing, a generalization has been introduced by Everett 
in [10]. Convergence and other properties have been further investigated in 
|15j |17j |20) |21j . As Renyi remarked in |2"T] , the formulation of Bissinger and 
Everett are somewhat stricter than necessary. In an elementary formula- 
tion (additional requirements could be imposed to obtain convergence), the 
so-called " f-expansion" of Bissinger and Everett can be translated into an 
expansion system as follows: 

Example 3.27 (f-Expansion). Let ^ be such that each [Si^Ui] equals 
([0, 1),0), and let / : [0, 1) — t- / be a monotonic mapping from [0,1) onto 
a certain interval I of the extended real line M U {— oo, +00} . Let each Ci 
equal Z U {—00, +00}. We set all Pi identical to 

P^y=[f{y)\. 

and all Ei identical to 

Eiy = f{y)-Py = f{y)-[f{y)\, 

where [— ooj := —00 and [+00J := +00. It should be noted that, because 
[0, 1) is an half open interval, the value +00 can only be attained for strictly 
decreasing /, and —00 for strictly increasing /. This expansion system 
models Everett's /-expansion for strictly increasing /, and Bissinger's /- 
expansion for strictly decreasing /. 
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4 Isomorphisms between Expansion Systems 



In this section we will introduce the concept of homomorphisms, and more 
importantly isomorphisms between ES's. When two ES's are isomorphic, 
then they essentially have the same underlying structure, and so they also 
share the other properties that depend on this structure. The isomorphism 
itself serves as a translation from one ES to another. Definition |4. 1| will give 
us the tools to construct such homomorphisms and isomorphisms. 



4.1 Homomorphisms and Isomorphisms 

Definition 4.1. Let ^' = {S,C,F) and = {S',C',F') be two expansion 
systems. A pair of maps A = (A'^,A'^) = ({Af {A^}jgN)i such that 



Af: Si 
A?: Ci 



C' 



(4.1) 
(4.2) 



is called a homomorphism of expansion systems, if the following equations 
hold: 



A9Pi 



P'A^ 



(4.3) 
(4.4) 
(4.5) 



If this A consists of only injective maps, then A is called an embedding of 
^' in ^f'. If A is bijective, then we call it an isomorphism and we say that ^ 
and ^I'' are isomorphic. 

Based on this definition, we can already get an idea of its purpose by the 
following proposition. 

Proposition 4.2. Suppose ^ and ^' are ES's, and A is an homomorphism 
from ^ into If y £ then A^y e S'J"^ and we have (A^y) = 



Proof. If y G , it follows by Theorem 



2.3 



that there is an element yl*'"! g 



Sq with the same coefficients as y up to and including order n — 1. Now by 



(|43l) and (0|, we have that ^„(A^ ?/[*'' 



and thus A^ yl*'*^! G S^ 



By (4.5) we have that Agy has the same coefficients as A^y up to and 
including order n — 1, and thus Theorem 2.3 tells us that A^y £ S'q"^ and 
(A^ y) =A^y[*'"l. 



□ 



The equations (4.4) and (4.5) can be rephrased in terms of diagrams. 



These requirements are satisfied, iff the diagrams in Figure [2] commute. 
From these diagrams it easily follows that if A is bijective, then the set of 
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5,: 



A? 



S' 



Si+i 



Si > Ci 



A? 



A? 



i+l 



S' 



Figure 2: Equations (4.4) and (4.5). 



inverse mappings A ^ = ((A'^) ^, (A"-") ^) is an isomorphism as well. This 



gives us the following corollary to Proposition 4.2 



Corollary 4.3. // ^ and ^' are isomorphic by the isomorphism A, then 
y G Sq^^ iff K^y G Sq "'\ and we have: 



y 



,[*',n] 



and 



(4.6) 
(4.7) 



Remark 4.4. In particular, this corollary says that if ^ and ^' are isomor- 
phic, then properness of one of them also guarantees properness of the other. 
This cannot be said of accurateness, as the concept of isomorphy is indiffer- 
ent to which topology is given on the initial element space Sq. However, if 
we have that Aq is a homeomorphism between Sq and Sq, then of course con- 
vergence or accurateness of one of the ES's, implies respectively convergence 
and accurateness of the other. 

Furthermore, it should be noted that the relationship of being isomorphic 
defines an equivalence relation among expansion systems, as is easily verified. 

We will now regard some elementary properties for commutative dia- 
grams, which we will then apply to the situation above. These properties 
can be best illustrated with the diagrams of Figure [3j The first of these four 



X 



Y' X 



Y' X 



Y' X 



A, 



Y > Z Y 

exists 



Z Y > Z 

unique 



Y' 



Y ' > Z 

(a) (b) (c) (d) 

Figure 3: Four Properties of Commutative Diagrams. 



diagrams, Diagram (a) , says that if we have only three of the four maps, and 
there is an injective map i leading to the map that is missing, then a fourth 
map (pa niust exist, which completes the commutative diagram. The same 
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can be said if there is a surjective map s, such that the map that is missing 



would lead to the map s, as is shown in Diagram |(b)[ Diagram (c) expresses 
that if we have three out of four maps, and a surjective map s leads to the 
map missing, then any map 0c completing the diagram has to be unique, 
i.e. there is at most one map (j)c which could complete the diagram. Finally, 
Diagram | (d) | says the same thing for the case that the missing map (/>rf would 
lead to some injective map i. One may verify that maps are a solution for 
these diagrams, if and only if we have: 

(l^aliix)- Kx)^ fix), (4.8) 

Mx) e s-Hf{x)), (4.9) 

(j)c: s{x)^ fix), (4.10) 

(pd- X ^ i,^^ o f{x), (4.11) 

where the map / : X — t- Z denotes the composed map from X to Z via 
Y' , and i^^ is the left inverse of i on f{X). The map <f>d in (4.11) can be 



defined if and only if f{X) C i(Y). For our situation of expansion systems, 
the above statements have, among others, the following implications: 

Theorem 4.5. Let us be given two sequences of element spaces S and S' , 
and two sequences of coefficient spaces C and C . Throughout the follow- 
ing statements, we do not make any prior assumptions about existence of 
the connecting sequences E,E' ,P,P' or IsP , except for those which are 
mentioned explicitly. Then we have that: 



Al Given E and injective A"^, then an E' satisfying (4.4) exists. 



A2 Given P ,ISp and injective , then a P' satisfying (4.5) exists. 



B Given E' and surjective IsP , then an E satisfying (4.4) exists 



CI Given E and surjective , then the unique map E[ satisfying (4.4) 
exists, iff it sends Af y to Af_^^£'j y. 



C2 Given P,hP and surjective , then the unique map P' satisfying (4.5) 
exists, iff it sends Af y to A^'Pi y. 



D Given E' and injective A"^, then the unique map E satisfying (4.4) 
exists, iffE[Kl{Si) C Af,,(5,+i). 



Remark 4.6. Different from the other statements of Theorem 4.5 , the claims 
B and D can be used in case the ES ^' is given, and one wants to embed a 
certain = {S,C,F) into this ^' . w 

Theorem 4.7. Let ^ be an ES, and let S' and C be two other sequences 
of respectively element and coefficient spaces. Then there exists an ES ^' = 
(S' ,C' , F') isomorphic to iff there is a collection of bijective mappings 
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between C and C , and bijective mappings between S and S' , satisfying 
(4.3). Given such bijective and JSp , F' is uniquely determined by: 



El = Al,EaA',y' 

pl = A^paA^y'. 



and 



(4.12) 
(4.13) 



Proof. The existence of a sequence F' = {P',E') satisfying (4.4) and ( |4.5[ ) 
fohows from Theorem 4.5 f\l and A2. As A'^ is bijective, Theorem 4.5pi and 



C2 tells us that E' and P' are uniquely determined by (4.12) and (4.13). □ 

Corollary 4.8. Let ^ = {S, C, F) be an ES. Then an ES = {S, C, F') 

is isomorphic to iff there is a collection of automorphisms AP on C, and 
automorphisms A^ on S, satisfying (4.3), such that F' can be defined as in 
( |4T2| ) and ^A^. 

Theorem 4.9. Let ^ be an ES, and suppose each map Ei is a composition 



of two maps E^''' : Si — ?■ S'^ and E^' : S[ — )• S'i+i, E\^' being bijective. Let 



^' = {S', C, F') be a another ES, where S[ is the domain of E\ , and which 
is further defined as: 



C[ 

Pi 
E' 



In diagram, this looks like Figure 
isomorphism from ^ to ^' . 



Ei Ui 

P^o{E^)-\ 
E^^oFf^. 



(4.14) 
(4.15) 

(4.16) 

(4.17) 



Then A^ = E^ and A^ = Id form an 
Proof. We observe that our A'^ and A*-^ are bijective, satisfying (4.3), and 



thus the claim follows from Theorem 14.71 



□ 



4.2 Applications of Isomorphisms 

We will now continue with some examples of isomorphisms between expan- 
sion systems. 

Example 4.10 (Newton's interpolation formula). Consider an ES ^' 
where each Si is the space of real polynomials, which we will view as func- 
tions on M in the variable and we let = 0. Let Pi be given by 
Pi{y) = y(0), and let Ei := A, which is the forward difference operator: 
Ay(x) = y{x -|- 1) — y{x). Now Fi = {Pi,Ei) is injective, for suppose we 
would have y(0) = z{0) and A(y) = A{z) for certain polynomials y,z G Si. 
Then by linearity A{y — z) = and (y — z){0) = 0, which is only satisfied by 
the zero polynomial. Hence y = z, and we conclude that ^ is proper. As we 
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Figure 4: Expansion systems with a relative shift 



have that /S!^y = iff /c is bigger than the polynomial degree of y, Definition 



2.1 tells us that the order of an element y equals its degree plus one (except 
for the order of 0, which is just 0). So all elements in Sq are elements of 
finite order, and thus by Remark |3.2| we have that ^ is an accurate ES, 
regardless of the topology on 5*0 . 

However, this theoretic result doesn't give us an idea of how these kind 
of approximations look like. To find this out, let's not write out our poly- 
nomials in the more common basis 1 , etc., but in the base 



X x{x — 1) x{x — l){x — 2) 
'1!' 2\ ' 3^ ' 



where the (^) denotes a binomial coefficient. Now if for certain i, y\^i 
happens to be of the form 



= Ci+i ... + Cn-i 

where Cj = Pi y, then using that 



X 



n-i-2 
Ci), we get: 



Fi [ci,yl^^) = Ci + Ci+i\\ + ... + Cn-i\ . ^ 



\ti\ 

That y\j^i is indeed expressible in the form above, follows by induction, 

[nl 

whereby we observe that yn fits in this pattern as being the empty sum 0. 
The n-th order approximation thus equals: 



Co + ci 



+ . . . + Cn-l 



n-1 
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As we had already proven that ^ is convergent, it follows that every poly- 
nomial y can be written as: 

k>0 ^ ' k>0 

This result is known as Newton's forward difference interpolation formula 

(cf. mi §1.2]). 

No isomorphisms so far. Now let's construct a new ES which element 
and coefficient spaces are exact copies of S and C. The maps 

Af y(x) = -y(-x), and (4.18) 
Afc=-c, (4.19) 

are then clearly automorphisms on Si and Cj (they are even self-inverse 
functions), and A'^ satisfies (4.3). Therefore, Corollary 4.8 tells us that if we 
put: 

Pi y{x) = Af (Af )-^ y{x) = Af - y{-x) = 

A? -y(0)=y(0), 

and furthermore 

Elyix) = AtiP.(Af)-iy(x) = Af^^P, - y(-x) = 

Af+i(-y(-x - 1) + y{-x)) = y{x) - y{x - 1), 

then ^' = (S, C, F') will be isomorphic with ^. Note that we here have 
P' = P, and E'- equals the backward difference operator V. So the only 
difference with is that our expansions are given by backward differences, 
instead of forward differences. To establish properness and accurateness of 
^' , and to find out what approximations it provides, although we could, we 
do not have to repeat all the steps we have gone through for Simply by 
invoking Corollary 4.3 and Remark 4.4, we conclude that ^' must be proper. 
Giving both and S'q the discrete topology, Aq becomes a homeomorphism. 
As we have already remarked, ^ is accurate for every topology on Sq, in 
particular the discrete topology. Therefore, has to be accurate as well, 
and for its approximations we find: 

AS( E (-^)*7V°) ....(.-..^1))^ ,4.0) 

fc=0,...,n-l 

(_l)fe+lV'=y(0) , , , , , 

^ kl {-^)---{-X-k + l) = 

fc=0,...,n— 1 

y: (4.21) 

k=0,...,n—l 
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where in (4.20) we use both that A{—y) 



■(Vy)(— x), which imphes by induction: A 



: —Ay and that A(y(— x)) = 
-y{-x)) = (-l)'=+i(V'=y)(-x). 



Taking the hmit of (4.21) as n goes to infinity, we get: 



V''l/(0) 



fc>0 



kl 



(x + k-l), 



(4.22) 



which is known as Newton's backward difference interpolation formula (cf. 
[m Ch. 1, prob. 6]). 

Exa mple 4.11. The expansions of the decimal expansion, described in Ex- 



ample 



1.5 



are composed of two maps E^^'^ and E^'^\ namely: 
eI''>: [0,1) ^ [0,10), y^lOy, 

E^^^: [0,10) ^[0,1), y'^y'-ly'\. 

As the first of these two maps is a bijection, we have by Theorem |4.9| that the 
is isomorphic with the ES described by: Si = [0, 10), C'^ = 
10(y' — -P/y). The elements and its convergents 



1.5 



^' of Example 
Ci,Piy'=[y'\ and Ely' 
from the former ES can be translated into elements and convergents for this 
ES, just by multiplying with a factor 10. 



Something similar applies to the continued fraction (Example 1.6), and 
we have that the ES is isomorphic to one defined by: 5- = M>i, C'^ = Ci, 
Piy'=[y'\ andE'y' = l/{y'-Ply). 

There are quite some more situations, for which there exists the choice in 
what order to apply the mappings which the expansions Ei our composed of. 



In those cases. Theorem 4.9 may tell us that the chosen order is actually not 
essential for the approximations which the ES provides. Another example 
of this can be found in the appendix about approximation systems. 
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Appendix: Approximation Systems 



In this appendix, we will have a glimpse at so-called approximation systems. 
A first description of approximation systems was given in |18j . and in more 
detail it is discussed in [19]. Here, we will give another description of AS's 
(approximation systems), but now in the language of expansion systems. 
The research regarding approximation systems is still going on, and there are 
quite some issues that need yet to be clarified. The content of this appendix 
should therefore be read as a rough account of recent investigations, rather 
than a rigorous treatment of the subject. 

Although the original formulation of AS's is somewhat different from the 
language of ES's, we will see that ES's are indeed capable of describing 
AS's, and even in a more general setting. However, this translation from 
AS's to ES's is not so straightforward, and particular complexities do arise. 
The different nature of AS's and ES's is the main reason for this. While 
the first is described by a sequence of differential equations for a particular 
sequence of functions, the second one has to be described by a sequence of 
maps between function spaces as a whole. 

Perhaps we could avoid some problems by drastically limiting the extent 
of the spaces Si, but besides that this would still leave some obstacles to 
overcome, it would also be a pity, as it turns out that the use of expansion 
systems allows us to give a general approximation algorithm, which can be 
applied to a wide range of function at the same time. The formulation of 
AS's chosen in this section is based on function germs. This solves certain 
technical problems, that could otherwise prevent the ES from being well- 
defined. But presumably, it is not strictly necessary to rely on the concept 
of germs, and at the end of this section we will briefly discuss the possibility 
of describing the element spaces of an approximation system as consisting 
of analytic functions on a certain domain. 

Instead of giving a complete definition of AS's, we will start give some 
examples of it. These examples will should already give an idea of what 
AS's are. 

Examples of Approximation Systems 

Before we come to the first examples, we mention some notation regarding 
function germs, which will be used in the examples. For a brief description 
of germs, see for instance |14l end of Ch. 3]. 

Definition 4.12. Given an analytic germ y at xq, we denote the multiplicity 
of its zero in xq by M{y). This can attain all values in Z>o, as well as oo, 
which we reserve for the germ constant 0. 

Example 4.13. Let Si be the space of germs in xq with constant term 1, 
and let Vi = 1. Further, we let Cj = C x (Z>o U {oo}), and let Pi be given 
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by 



Piy 



-M{Dy) 



Dy 



(xo) , M{Dy) 



where the first entry simply denotes the first Taylor coefficient after the 
constant term 1. If y is constant 1, we may set Piy = (0, oo) and Eiy = \ 
(satisfying (1.1)). For all other y we set: 

Dy{x) 



Ei{y) 



Ci{x - XoY 



(4.23) 



where (cj, rrii) = Pi y, and Oi can be any value in C\{0} (note that the germ 
h := Dy{x) / {ci{x — xq)"^^) is 1 in xq, so that we can take the principle value 
for the logarithm in = exp(ajlog h) ). These mappings define a proper 
ES ^, for we have that Fi = {Pi, Ei) is bijective, with inverse: 



Ff' (c, m,y) = l+ c{t- xoTy{t) dt, 

Jxo 

which evaluates to 1 in case c = 0. In the context of our ES, all these 
expressions are formally to be interpreted as just germs in xq, but as we will 
see later, they have more practical value than just representing germs. 

Remark 4.14. The expansions Ei can be decomposed into E^^^ : Si — )• S'^ and 
eP: S'^S^+^, where: 



Ef'y' 



Dy, (4.24) 

^c,(x-xo)™»y ^ ' 

As the first of these maps is bijective on the space of germs with constant 
term 1, we have by Theorem 4.9 that our ^ is isomorphic with where: 



E'i y{x) 



( (x-xo) 



-Miy) 



y 

y{x) 



[xo) , M{y) 



[{x - xo)*^(^)y] (xo) 



Note that we have: 



p, = p;eP. 



(4.26) 
(4.27) 

(4.28) 



In fact, the formulation of ^' has been the underlying principle in |T8] . 
whereas the formulation of ^ coincides with the formulation in [TH]. As 
y = Dy gives an isomorphism between the spaces Si and 5-, and S'^ 
does not have the restriction of having constant term 1, it is clear that this 
restriction on Si is actually not that limitative. In fact we could drop this 
restriction on Sq (i.e. we take as initial space S'q), and in a similar way as 
we have mentioned in Example |1.5[ we let our first expansion be given by 

(2) 

simply Eq , keeping all further expansions just as they were. Then the ES 
can be applied to any holomorphic germs. 



35 



As we haven't fixed the value of in the AS in Example 4.13, it actually 
represents a whole collection of AS's. We will now consider some particular 
instances of it. When each = 1, it can easily be verified that the approx- 
imations are just the Taylor approximations (besides the restriction in 
on the constant term, the main difference with the approximations provided 



by Example 2.6, is that here the n-th order convergent gives the first n — 1 



non-zero terms of the Taylor series, whereas in Example 2.6 the terms up to 
and including order n — 1 will be given) . 

More generally, when each is a positive unitary fraction 1/pi with 
Pi £ Z>i, then the resulting approximations will be polynomials. Using 



concise notation, such as we have used in (1.7) and (3.9), we can write the 
approximation process as: 

y{w) = I + CQWi""' (^1 + ' ciW2'^' . y dw2^ dwu (4.29) 

using local coordinates w = x — xq, and where {mi,Ci) = PiU- This ex- 
pression suggest that this ES is actually accurate. Of course, by Theorem 



3.6 we have convergence in the weakly discrete topology. In fact, it can be 
shown that this means that the first n — 1 non-zero Taylor coefficients of 
the n-th convergent coincide with the first n — 1 non-zero Taylor coefficients 
of the initial function y (cf. |19| Cr. 5.4]). This means we have also weak 
convergence in the sense that the derivatives of all orders convergence in the 
point Xq. But in practice we often have even more than that. Consider the 
situation in which each pi = 2, and look at the element y{x) = 1/x and 



Xq = 1. Then (|4.29|) reduces to: 
1 




^. . dx^j (ix2 j dxi. 



(4.30) 

The coefficients Cj in this formula, have been moved outside of the integrals. 
The first few convergents are given by: 





= 1 








yW(x) 


= 1 + 


X 

2 






yl'kx) 


= 1 + 


X 

2 + 


3x2 

~8~ 


3x3 


y^'Hx) 


= 1 + 


X 

2 + 


3x2 
8 


5x3 



175x4 147x5 343x5 343x5 



1024 2048 16384 131072 

Interpreting these polynomials as complex functions rather than germs, con- 
vergence can be established on the open disk with radius 1. In [19] this has 
been proven for all powers p > 2. However, despite the fact that the head 
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of all these polynomials coincides with the Taylor series, computer simu- 
lations have pointed out that the actual domain of convergence is in fact 
larger than this open disk, and thus not symmetric. Of course, the domain 
of convergence is symmetric in the complex direction, and the pole in x = 1 
lies on its boundary. What the exact contours are, remains nonetheless a 
challenging question. 

Remark 4.15. In the example above, if Si would be defined as a space of 
holomorphic functions on a fixed neighborhood U of xq, then we can always 
find a function y on which Ei is not well-defined, as it would introduce a 
branch point to U. And at the same time, possibly this y could be success- 
fully approximated, if only our domain was somewhat smaller than U. In 
the formulation with germs, the algorithm applies without problems to all 
relevant functions. In the subsection Further Aspects of Approximation Sys- 
tems we will briefly discuss an alternative formulation, which could possibly 
solve this issue in a yet different way. 



We have allowed the ai in Example 4.13 to depend on i, and it is indeed 



possible to vary as i increases. For example, when we let a, run through 
2,5,2,... and we look at y{x) = e^, then we get the expansion: 

e^ = l + £(^l + i£^(l + i£^(...)V)^x,) dx, (4.31) 

Although the convergence of this formula turns out to be very fast, it also 
requires much more effort to evaluate than a Taylor series. 

So far we have only considered unit fractions, which give polynomial 
convergents. Another possibility are positive integer a^. For instance, when 
we let ai run through 2, 3, 4, . . . instead, and pick y = again, then we get 



more or less the opposite of (4.31): 



= 1 + 




1-^2!/ {/ 1 + 3! / (/. . . dx3 dx2 dxi. (4.32) 



Here, the convergents are no longer polynomials, but more complex expres- 
sions. Symbolic computation is impossible for higher orders of the conver- 
gents, and so the only option seems to be numerical evaluation, or computa- 
tions in terms of power series. Although the appearance of the factorials in 
the last two formulas for may seem remarkable, they are actually imposed 



by our choice of Oj. For example, if we would keep ai constant as in (4.30) 
then the coefficients form a geometric sequence instead. 



Another interesting instance of Example |4.13 we get by setting ai = —1 
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In concise notation, this gives us expansions like: 

1 

1+ " 



1 



(4.33) 







1 + 



X2 



dX3 



X 



1 + 



1 - 

a dxi 



X3 



1 + 



1 — a dx2 



(4.34) 



1 + 



X2 



a dx-t 



cosh X 



1 + 



f 

Jo 



1 + 

Xi dxi 



'^'3 1 — a dx4 



Xl 



X2 dx2 



(4.35) 



3 + 



X2 



7 X3 dX3 



^•3 221X4^X4 

2205 • . . . 



Remark 4.16. The first two of these formulas show a regular pattern. In fact, 
we can view these two formulas as a repeated application of an operator T, 
for which y is a fixed point. For instance, in (4.33) we have: 

Ty{x) 



1 + 



1 



dxi 

'■^'^ dx2 

yix2) 



which has as fixed point . The reason that there exists such an operator T, 
is that the initial function y satisfies a finite cycle of differential equations, 
determined by (4.23). In this case, it is a cycle of length 2. In general, if we 
are given a cycle of equations: 

yiix) = Ci+ fi{yi+i{t),t)dt, (i = 0, . . . , n - 1, y„ = yo), (4.36) 

then the solution may also be found by a Picard iteration: 

Ty = c+ [ foy{t)dt, 

on a multi-dimensional y : C — > (cf. [T^", 7.2]). For this, we set / : ^ 
as: 

fiVO, ■ ■ • ,yn-l) = (/O o yi, • • • , fn-2 O yn-1, fn-1 ° Vo) 

and as constant of integration simply c = (cq, ci, . . . , c„_i). 
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Starting this iteration with the vector y = {vq, . . . (for the case 

each = 1, this is just the vector (1, . . . , 1)), then after n repetitions, the 
first entry of this vector equals y'"!. 

Of course, the algorithm still applies when there is no such cycle for y. 
The third formula is an example of this, and its coefficients do not follow an 
obvious pattern. Such continued fraction-like expressions can be obtained 
for every analytic function, and for most it will require quite some effort to 
calculate the coefficients. Furthermore, the convergents cannot be expressed 



in terms of elementary functions for high orders of n. For example, in (4.33) 
we can only determine the convergents up to and including order 4: 

y[°l(x) = 1 

yW(x) = 1+ X 

y[2l(a:) = 1 -log(l - j;) 

y[3l(x) = l + e-h^ - e-h— , 

e e 

where li denotes the logarithmic integral. For the orders n > 4, we have 
to rely on numerical methods to evaluate the expression. Although the 
rate of convergence of the above formulas turns out to be rather slow, its 



domain of convergence is surprisingly large. In the situation of (4.34), the 
approximations even seem to converge when we evaluate them on paths 
going around the branch point of x", when a is not a nonnegative integer. 
We will come back to this observation at the end of this appendix. Before 
that, let us consider two more examples of AS's. 



Example 4.17. Unlike Example 4.13[ we now let Si be the space of holo- 



morphic germs at xq with constant term 0, and = 0. We let and Pi be 
just as in Example |4. 13 , but now we set y = if y = 0, and otherwise: 



These mappings define a proper ES ^, for we have that Fi = {Pi,Ei) is 
bijective, with inverse: 

Ff^c, m, y) = r c{x - xo)'"e2'W dt, 

J XQ 

which evaluates to in case c = 0. 
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This AS provides us with formulas such as: 



log X 



tan X 



i: 



dxi e 



^1 -I dxs e 
dx2 e ■'1 



dxA e " 



X I dx2 2x2 e 

dxi eJo 



f 

Jo 



dxs 0X3 eJo 



dxi j|x4 c 



(4.37) 



(4.38) 



For clarity, we put the variable of integration in front of the integrand, 
instead of behind. The first of these formulas can be viewed as iterating 



a single operator, like we saw with (4.33) and (4.34). In this case, the 



operator entails only one integration in total, as the cycle of differential 
equations consists merely of the equation: 



dx 



logx 



- log X 



In other words, this formula just expresses a Picard iteration. Note that in 
Equation 4.37, we have multiplied both sides of the equation by a factor — 1 



afterwards. The second formula does not arise from a cycle of differential 
equations, and so the coefficient do not follow a cycle either. Numerical 
simulations suggest that both of these formulas converge on a certain domain 



(as for (4.37), which is equivalent to a Picard iteration, this should not come 
as a surprise) . 

Example 4.18. In this example, we let Si consist of the germs with constant 

but now we let Ci = X (Z>o U {00}). 



4.13 



term 1 again, as in Example 
Further, we let Ei be the composition of maps E^-^'* and Ef'\ where E"-^^ 



is defined as in (|4.25|), and E'-^^ is the composition K o D, where Ky :- 



y — y{xo). This basically means that E^^'^ first differentiates a function, and 



chops off the constant term afterwards. Note that this E, 
any longer. We set Pi to: 



(1) 



isn't bijective 



Piy = {Dy{xo), (x'^^^ ^^^i^ Dy) (xo) , M{K Dy) ) . 

These mappings define a proper ES and Fi = {Pi,Ei) is bijective, with 
inverse: ^ 

Fi-^{b,c,m,y) = l + bx+ / c{t - xo)'^y{t)^^ dt. 



We may verify that this expression is mapped to y by Ei, and that the 
projection Pi indeed gives the parameters {b,c,m). 



The above AS gives us among others the expansion: 
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+ = 




^ . . . dx^ dx2 dxi. 



(4.39) 



This equation is obtained by applying the AS to (x + l)^, and subtracting 
the equal terms (i.e. 1 + 3x) afterwards. The coefficients obey the pattern 
bi = 3/2' and a = 3/2'^'''^. 

Further Aspects of Approximation Systems 

We have just seen some examples of approximation systems, as well as the 
kind of results it gives. This may have sparked a couple of questions. First 
of all, what is the meaning of all these formula's? As noted before, when we 
interpret the equations as dealing with germs, the equality symbol expresses 
that the left hand side is the limit of a sequence defined by the right hand 
side, namely in the weakly discrete topology. But in many cases, the expres- 
sions do also converge in the more common topology. For some particular 
cases, this has been established in [19j . For many others, this has only been 
suggested by numerical simulations. 

Although for some AS's the rate of convergence can be really high with 
respect to particular elements y, in most cases this does not seem to com- 
pensate for the effort necessary to (numerically) evaluate these kind of ex- 
pansions. Of course, it may well be that there are situations in which the 
algorithm really serves practical purpose, but in other situations it would 
probably not be the prime motivation to study these identities. Probably 
more interesting is the fact that these expansion do converge in the first 



place, also when it involves quickly increasing functions such as in (4.38). 



Even more remarkable, is that the domain of convergence can take extraor- 



dinary shapes. We already mentioned this regarding (4.30) and (4.34). 

For the first of these, we know for sure that its domain of convergence is 
at least the open disk with radius 1 (see [I9j), but evaluating the algorithm 
numerically, we find a domain which is considerably larger than that, delim- 
ited by a very irregular boundary, with even fractal like characteristics. As 
for now, we could only speculate about what this is supposed to mean. 



At Equation (4.34), we already mentioned that here the approximations 
seem to converge, even when going around the singularity at x = for 
a ^ Zi>o. What we mean with this is the following. Let's first choose 
a parametrization 7: [0, r] — t- C \ {0} of a path going around the origin 
(this could be ^{t) = e** for example, but if we are performing a numerical 
simulation, a path in the shape of a square would also do). Next, we evaluate 



the expansion given in (4.34) up to a certain order, but when calculating the 



integrals, we restrict ourselves to the path 7. Then, we determine the values 
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of x°- on the same path 7 by analytic continuation (if 7(t) = e*% this is just 
e"*'). If our 7 doesn't happen to have crossed a singularity of some (of 
which there can only be finitely many), it turns out that the approximated 
values of y^"^ on 7, converges to the value of y on this 7, as n gets larger. In 
somewhat loose notation, this means: 



^0 



a D'y{ti) dti 



1 + 



/' 

Jo 



1 - a D-fit2) dt2 



1 + 



t2 



a D^{h) dts 



;*3 l-aDj{U)dt4: 



/ 

Jo 



where the integration variables U run through the domain [0, r]. As long 

as 7 doesn't accidentally cross some singularity like we just mentioned, this 
identity always seems to hold, no matter what the end point of 7 is, and 
how many times it has turned around the singularity at x = 0. Performing 
a similar procedure on the Taylor expansion (viewed as an AS), we don't 
have this property. For example, picking ^(t) = e*', we can find either by 
calculation or by reasoning, that: 

7(t)« = 1 + a(e^* - 1) + ^^^^(e^* - 1)' + . . . 

which diverges for i > |. 

In order to study this kind of convergence behavior for AS's in general, 
it may be desirable to give a different formulation of AS's. Possibly, we 
could let the spaces Si consist of multilayered functions covering the com- 
plex plane, allowing different functions to live on different domains, but each 
of these domains being marked with a certain basepoint xq. Then the expan- 
sions Ei determine a transition from one multilayered function to another, 
which might as well alter the Ricmann surface a function lives on. If we 
want to compare the values of two functions y and y' that live on different 
domains, (so that we are able to define a criterion for convergence), we can 
lift paths in C, such that the begin point of the lift is in both situations xq- 
The difference in value of y and y' on the lift of a suitable path 7, serves as 
a measure for how close the two functions are on that particular path. Of 
course, to give this whole formulation a solid foundation, will require a lot 
more work, which we leave outside the scope of this thesis. 

After all this discussion, there is at least one question which still isn't an- 
swered yet: what is an approximation system? The different examples may 
of course have given an impression of what is meant with it, but a strict 
definition has not been given. As one may have noticed, the expansions and 
projections take various forms, which are difficult to capture all in a single 
definition. Restricting ourselves to the formulation with germs, approxima- 
tion systems may be informally described as expansion systems on spaces 
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of holomorphic germs, for which each expansion is a transformation which 
(possibly) erases some of the information of the lower order derivatives (i.e. 
an initial part of the Taylor series), and for which the projections are the 
containers of this information. A further smoothness condition on the ex- 
pansions and projections may be desirable, in order to have convergence of 
the AS. Of course, this description of AS's is still to vague to be used within 
a mathematical context, but hopefully it does give an idea of the mechanism 
behind AS's. 

It is clear that approximation systems open the door to a virtually endless 
number of expansions. And then we haven't even mentioned yet the case in 
which we do not involve any differentiations in our maps Ei. For instance, 
if we replace the differential operator in Example |4.13 and Example 4.17 by 



the operation K y = y — y{xo), as introduced in Example 4.18, then a whole 
new family of identities arises, such as: 

W{x) = xe~^^ (4.40) 



2 + x 



^ 1 + x-^ 




(4.41) 



tan X 



X 



(4.42) 

where W{x) denotes the Lambert W function. In the literature, the first 
two kinds of expansions are respectively referred to as infinite exponentials 
(cf. [2] ) , and nested radicals (cf . [1] ) . The third type of expansion is again a 
continued fraction, but different from Example |1.6[ now including a variable 
x. For arbitrary analytic functions, the approximation provided by this ES 
corresponds with the so-called C-fraction for a given analytic function, or 
even formal power series (cf. [T^l 2.4]). This particular identity in (4.42) 
belongs to the class of Gauss 's continued fractions, and is more specifically 
called Lambert's continued fraction (cf. [23t 91.7]). 
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